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METHOD OF DETECTING SPATIALLY CORRELATED VARIATIONS 



IN A PARAMETER OF AN INTEGRATED CIRCUIT DIE 



FIELD OF THE INVENTION 



The present invention relates generally to the 



testing of integrated circuit dies on a wafer during 
manufacture. More specifically, but without limitation 
thereto " the 'plresentT invention relates to reducing the 
variation of a selected parameter in a production lot of 
integrated circuit die. 



integrated circuits is detecting and rejecting integrated 
circuit die replicated on a silicon wafer that exhibit 
values of a selected parameter, for example, quiescent 
current (Iddq) , that differ significantly from a mean value 
of the parameter. Integrated circuit die having values of 
the selected parameter that differ from the mean value by 
more than a selected threshold are called statistical 
outliers. Statistical outliers may pass performance 
testing, however, they may be more subject to premature 
failure and thus reduce the average service life of a 
production lot. Accordingly , a need exists for a method 
of detecting the statistical outliers. 



BACKGROUND OF THE INVENTION 



An important issue in the manufacture of 
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SUMMARY OF THE INVENTION 



In one aspect of the present invention, a method 



of detecting spatially correlated variations includes 
5 measuring a selected parameter of each of a plurality of 
electronic circuits replicated on a common substrate; 
calculating a difference between a value of the selected 
parameter at a target location and that of an identical 
relative location with respect to the target location for 

10 each of the plurality of electronic circuits to generate a 
distribution of differences; calculating an absolute value 
of the distribution of differences; and calculating an 
average of the absolute value of the distribution of 
differences to generate a residual for the identical 

15 relative location. 



process for reducing the variation of a selected parameter 



parameter of each of a plurality of integrated circuit die 
20 replicated on a wafer substrate; calculating a difference 
between a value of the selected parameter at a target 
location and that of an identical relative location with 
respect to the target location for each of the plurality of 
integrated circuit die to generate a distribution of 
25 differences; calculating an absolute value of the 

distribution of differences; calculating an average of the 
absolute value of the distribution of differences to 
generate a residual for the identical relative location 
that is representative of an expected value range of the 
30 selected parameter at the identical relative location; and 
rejecting any of the plurality of integrated circuit die 
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In another aspect of the present invention, a 



of an integrated circuit die includes measuring a selected 



having a value of the selected parameter that lies outside 
the expected value range . 

DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

5 

The present invention is illustrated by way of 
example and not limitation in the accompanying figures, in 
which like references indicate similar elements throughout 
the several views of the drawings, and in which: 
10 FIG_. 1 illustrates a typical distribution of 

average die quiescent current for a replicated pattern of 
integrated circuit die on a silicon wafer; 

FIG. 2A illustrates a plot of residual quiescent 
current as a function of relative die location according to 
15 an embodiment of the present invention representative of a 
smooth process variation; 

FIG. 2B illustrates a three-dimensional mesh of 
the plot in FIG. 2A; 

FIG. 3A illustrates a plot of residual quiescent 
20 current as a function of relative die location according to 
an embodiment of the present invention representative of a 
non- continuous process variation; 

FIG. 3B illustrates a three-dimensional mesh of 
the plot of FIG. 3A; 
25 FIG. 4 illustrates an ordered rank of the top 4 0 

die that are the most highly correlated with the target die 
from the location average results of FIGS. 3A and 3B; 

FIG. 5 illustrates average die quiescent current 
plotted as a function of wafer location with 
30 stepper/reticle effects; 
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FIG. 6 illustrates a flow chart of a method of 
detecting variations in a spatially correlated parameter in 
accordance with an embodiment of the present invention; and 

FIG. 7 illustrates a flow chart of a process for 
5 reducing the variance of a selected parameter of an 

integrated circuit die replicated on a silicon wafer in 
accordance with an embodiment of the present invention. 

Elements in the figures are illustrated for 
simplicity and clarity and have not necessarily been drawn 
10 to scale. For example, the dimensions of some of the 
elements in the figures may be exaggerated relative to 
other elements to help to improve understanding of 
embodiments of the present invention. 

15 DESCRIPTION OF THE ILLUSTRATED EMBODIMENTS 

Current approaches for detecting statistical 
outliers are based on the assumption that die located 
closest to one another on the wafer have the most highly 

20 correlated parameter values, resulting in a dominant 

spatial pattern that is continuous. In other words, the 
closer the die are to each other on the silicon wafer, the 
more highly correlated are the corresponding values of a 
selected parameter. However, the dominant spatial pattern 

25 may not be continuous. For example, die located at similar 
reticle positions may result in corresponding values of a 
selected parameter that are less correlated for neighboring 
die than for die located at similar reticle positions. 

A reticle is a finely resolved picture that often 

30 contains multiple images of the same die pattern. The die 
patterns in the reticle are transferred simultaneously from 
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the reticle to the silicon wafer according to well known 
techniques of photolithography as the reticle image is 
stepped across the silicon wafer. Due to the effects on 
the light traveling through the reticle to the silicon 
5 wafer, the exact pattern is not perfectly maintained across 
all copies of the die pattern. The result of these effects 
are often related to the relative positions of the die 
patterns within the reticle. If these effects dominate a 
process step that would otherwise cause a continuous change 

10 in a pairameter value, then the reticle dependence pattern 
results in a higher correlation of parameter values in non- 
adjacent die than in adjacent, or ''nearest neighbor", die. 

Using test data, die-to-die correlations of a 
selected parameter may be calculated by averaging measured 

15 values of a selected parameter of die having an identical 
relative location with respect to a target location of each 
die on the silicon wafer respectively. A target location 
is the location of a specific die on the silicon wafer. A 
relative location is the location of a die displaced from 

20 the target location by a specific difference in the X-Y 

coordinates. An example of a relative location is (-5,10), 
that is, 5 units in the negative-X direction and 10 units 
in the positive-Y direction with respect to the target 
location. For a target die having X-Y coordinates 

25 (150,300), the corresponding die having the relative 

location (-5,10) would have the X-Y coordinates (145, 310). 
For a target die having the X-Y coordinates (620,222), the 
corresponding die having the identical relative location 
would have the X-Y coordinates (615,232). For a target die 

30 having the X-Y coordinates (500,800), the corresponding die 
having the identical relative location would have the X-Y 
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coordinates (495,810), and iso on. For each target die on 
the wafer, there is a corresponding die having the 
identical relative location, except near the edges of the 
wafer, where the identical relative location may not lie on 
5 the wafer. 

Because a single production lot provides test 
data that allows for calculating die-to-die correlations of 
a selected parameter for each relative location from 5,000 
to 10,000 or more times, an adequate sample is provided to 
1 0 obtain a^ highly r^sqlve<^ of the 

selected parameter as a function of relative location. Die 
having values of the selected parameter that exceed an 
expected value range may be identified as statistical 
outliers and rejected from the production lot, thereby 
15 increasing the average service life of the production lot. 

By calculating die-to-die correlations of a 
g selected parameter of an integrated circuit die as a 

function of identical relative location on a silicon wafer, 
a highly resolved image may be obtained of not only the 
20 major spatially correlated patterns across the wafer, but 
also of more subtle patterns resulting from the 
manufacturing process. 

FIG. 1 illustrates a typical distribution of 
average die quiescent current for a replicated pattern of 
25 integrated circuit die on a silicon wafer. Shown in FIG. 1 
are target die locations 102 and corresponding die relative 
locations 104 . A gray scale is used to indicate the 
variation of the quiescent current at each of the target 
die locations 102. The variable gr is an index that 
30 identifies the specific relative location 104 of a 

corresponding die with respect to each of the target die 
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locations 102. The variable gr ranges from 1 to n, where n 
is selected to include as large or as small an area around 
each of the target die locations 102 as desired for 
correlating values of the selected parameter. In the 
5 illustrated example, n equals 49, and the selected 
parameter is quiescent current ( Iddq) , however, other 
values of n and other selected parameters may be used to 
suit specific applications , 

The residual of the selected parameter is 

10 calculated by measuring^ the value of the selected parameter 

at each of the target die locations 102 on the wafer and 
calculating the difference between the value of the 
selected parameter of a target die location 102 and the 
value of the selected parameter of the corresponding die at 

15 an identical relative location g with respect to each of 
the target die locations 102 to generate a distribution of 
differences. For example, if the measured values of the 
selected parameter at five of the target locations 102 are 
1.5, 1.2, 1.3, 1.6, and 0.9, and the measured values of the 

20 selected parameter at an identical relative location (-5, 
10) with respect to each of the five target locations 102 
are 1.3, 1.1, 1.4, 1.5, and 1.3, then the distribution of 
differences would be 1.5 - 1.3 = 0.2, 1.2 - 1.1 = 0.1, 1.3 
- 1.4 = -0.1, 1.6 - 1.5 = 0.1, and 0.9 - 1.3 = -0.4. The 

25 average of the absolute value of the distribution of 

differences for an identical relative location gr is defined 
as the residual of the relative location gr. The absolute 
value of the distribution of differences in this example is 
0.2, 0.1, 0.1, 0.1, and 0.4. The average of the absolute 

30 value of the distribution of the 5 differences in this 
example is 0.2 + 0.1 + 0.1 + 0.1 + 0.4 / 5 = 0.18. The 
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residual of the selected parameter for the relative 
location (-5, 10) is thus 0.18. In like manner, the 
residual of each of relative locations 104 may be 
calculated. 

5 FIG. 2A illustrates a plot of residual quiescent 

current as a function of relative location according to an 
embodiment of the present invention representative of a 
smooth process variation. Shown in FIG. 2 A are a target 
die location 102 and corresponding die locations 104 . 

^ 10 In FIG. 2A, a gray scale is used to indicate the 

O relative value of the residual quiescent current calculated 

H for each of the identical relative locations 104, where 

J white indicates the lowest variation in the quiescent 

current with respect to the target die location 102. As 
15 shown in FIG. 2A, the variation in residual quiescent 
current increases with distance from the target die X-Y 
location 102, which identifies smooth process variation as 
the only spatial pattern. In this case, a nearest neighbor 
residual process may be used according to well known 
20 techniques to identify statistical outliers, 

FIG. 2B illustrates a three-dimensional mesh of 
the plot in FIG. 2A. Shown in FIG. 2B are a target die X-Y 
location 102 and neighboring die locations 104. 

In FIG. 2B, the gray scale is further enhanced by 
25 adding a third dimension that emphasizes the spatial 
pattern of the residual quiescent current. 

FIG. 3A illustrates a plot of residual quiescent 
current as a function of relative die location according to 
an embodiment of the present invention representative of a 
30 non- continuous process variation. Shown in FIG. 3A are a 
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target die location 102 and corresponding die locations 
104 . 

In FIG. 3A, a gray scale is used to indicate the 
relative value of the residual quiescent current calculated 
5 for each of the identical relative locations 104, where 
white indicates the lowest variation in the quiescent 
current with respect to the target die location 102. As 
shown in FIG. 3A, the variation in residual quiescent 
current does not change monotonically with distance from 
10_ the target die locatJ^n 102 as in the example of FIG. 2A. 
The spatial pattern of the residual quiescent current in 
this example is non- continuous, therefore the nearest 
neighbor residual process, although possible, may be 
significantly improved for identifying statistical 

s 

15 outliers. A more accurate method of identifying 

statistical outliers is to use the mean value of the 
selected parameter as a function of relative location to 
determine a best estimate of the expected parameter for 
each relative location. 
20 FIG. 3B illustrates a three-dimensional mesh of 

the plot of FIG. 3A. The description of FIG. 2B also 
applies to FIG. 3B. 

Plots similar to those of FIGS. 2A, 2B, 3A, and 
3B may be generated to identify spatially correlated 
25 parameters in other electronic circuits replicated on a 

common substrate and in a variety of applications in which 
a common element is replicated in a spatial pattern. 

FIG. 4 illustrates an ordered rank of the top 4 0 
die that are the most highly correlated with the target die 
30 from the location average results of FIGS. 3A and 3B. 

Shown in FIG. 4 are a target die 402 and corresponding die 
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relative locations 4 04 numbered by correlation rank with 
the target die 4 02. 

As shown in FIG. 4, the relative locations most 
highly correlated with the target die lie on the same 
5 horizontal line as the target die. 

FIG. 5 illustrates average die quiescent current 
plotted as a function of wafer location with 
stepper/reticle effects. Shown in FIG. 5 are a target die 
502 and corresponding die relative locations 504. 

10 The _gray _scale_indicates the degree_ of 

correlation of average die quiescent current of the 
corresponding die relative locations 504 with the target 
die 502. As shown in FIG. 5, there is a combination of 
smooth variation and row dependent variation in quiescent 

15 current. 

FIG. 6 illustrates a flow chart 600 of a method 
of detecting variations in a spatially correlated parameter 
in accordance with an embodiment of the present invention. 

Step 602 is the entry point for the flow chart 

20 600- 

In step 604, a selected parameter of each of a 
plurality of electronic circuits replicated on a common 
substrate is measured according to well known techniques . 
The electronic circuit may be, for example, an integrated 

25 circuit die. The invention may also be practiced with any 
implementation of an electronic circuit or other device 
such as a micro-machine that may be replicated on a common 
substrate, and the term ''electronic circuit" as used herein 
includes all such types of devices including micro-machines 

30 that may be replicated on a common substrate. 
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In step 606, a difference is calculated between a 
value of the selected parameter at a target location and 
that of an identical relative location with respect to the 
target location for each of the plurality of electronic 
5 circuits to generate a distribution of differences. 

In step 608 , an absolute value of the 
distribution of differences is calculated. 

In step 610, an average of the absolute value of 
the distribution of differences is calculated to generate a 
10 representati^ye^ ya]^^^ for the res_idual _for the identical 
relative location. The calculated value is the mean 
absolute error, or if the median is used, the median 
absolute error. 

Step 612 is the exit point for the flow chart 

15 600. 

The method illustrated in FIG. 6 may be used to 
detect spatial correlation of a va'riety of parameters of 
practically any device that may be replicated on a common 
substrate . 

20 FIG. 7 illustrates a flow chart 700 of a process 

for reducing the variance of a selected parameter of an 
integrated circuit die replicated on a silicon wafer in 
accordance with an embodiment of the present invention. 

Step 702 is the entry point for the flow chart 

25 700. 

In step 704, a selected parameter of each of a 
plurality of integrated circuit die replicated on a wafer 
substrate is measured according to well known techniques. 
For example, the selected parameter may be quiescent 
30 current (Iddq) . 
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In step 706, a difference is calculated between a 
value of the selected parameter at a target location and a 
that of an identical relative location with respect to the 
target location for each of the plurality of integrated 
circuit die to generate a distribution of differences. 

In step 708, an absolute value of the 
distribution of differences is calculated. 

In step 710, an average of the absolute value of 
the distribution of differences is calculated to generate a 
representative value for the identical relative location 
having an expected value range of the selected parameter at 
the identical relative location. The expected value range 
may be, for example, the value of the selected parameter at 
the target location plus or minus one-half the residual. 

In step 712, integrated circuit die at identical 
relative locations having a value of the selected parameter 
that lies outside the expected value range are rejected 
from the production lot. In the example described for FIG. 
1, the die at the relative location (-5, 10) with respect 
to the first and the last of the five target locations 
would be rejected as statistical outliers, since they have 
values of the selected parameter that lie outside the range 
1.3 + 0.18. Other values for the expected value range may 
be used to suit specific applications. 

Step 714 is the exit point for the flow chart 

700 . 

The method illustrated in FIG. 7 may be used to 
reduce the variation in a selected parameter in a 
production lot of integrated circuit dies, thereby 
increasing the reliability of the production lot. 
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In another embodiment of the present invention, 
the lot averaging may be performed for each wafer X-Y 
coordinate so that a new set of best estimates is re- 
calculated for each X-Y position. Re- calculating the best 
5 estimate locations includes the effects of edge die versus 
center die on the wafer. 

In yet another embodiment, other wafers in the 
lot may be utilized for the location averaging. In this 
case, the best estimates for a given X-Y location may be 

10_, the identical location on another waf er in_the lot^ T^^^ 

technique may be improved by re -ordering the wafers in the 
sequence in which they were processed to ensure more 
accurate estimation. 

While the invention herein disclosed has been 

15 described by means of specific embodiments and applications 
thereof, other modifications, variations, and arrangements 
of the present invention may be made in accordance with the 
above teachings other than as specifically described to 
practice the invention within the spirit and scope defined 

20 by the following claims. 
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